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SUMMARY 


Protection  factors  were  calculated  for  cylindrical  concrete  barriers 
with  an  interior  partition  using  Monte  Carlo  procedures.  The  results  were 
compared  with  protection  factors  calculated  using  the  Engineering  Method, 
in  order  to  evaluate  the  ability  of  the  Engineering  Method  to  predict  pro¬ 
tection  factors  for  cylindrical  barriers  with  interior  partitions.  The 
Monte  Carlo  results  were  also  analyzed  to  give  a  better  understanding  of 
the  effects  of  the  interior  partition  and  to  aid  in  further  development  of 
the  Engineering  Method. 

The  cylindrical  barriers  were  of  concrete  and  were  infinite  in  height. 
The  radius  of  the  external  barrier  was  10.0  ft,  and  the  interior  partition 
radii  were  2.5,  5.0,  and  10.0  ft.  The  barrier  thicknesses  were  20,  40,  and 
80  psf  (lbs/ft2),  and  the  interior  partition  thicknesses  were  20  and  40  psf. 
The  source  was  an  infinite  plane  60Co  source,  assumed  to  represent  fallout. 

The  results  from  previous  Monte  Carlo  protection  factor  calculations 
for  20-,  40-,  and  80-psf  cylindrical  barriers  without  interior  partitions 
were  used  in  the  comparisons  to  present  barriers  with  interior  partitions 
of  zero  thickness.  The  barrier  configuration  and  source  were  the  same  as 
those  used  in  the  protection  factor  calculations  with  interior  partitions. 
The  previous  results  were  also  used  to  represent  barriers  with  a  pseudo 
partition  at  a  radius  of  10.0  ft. 


The  flux  energy  and  exposure  rate  angle  distributions  were  also  calcu¬ 
lated  at  a  height  of  3  ft  inside  the  barriers  with  the  interior  partitions. 
These  data  were  obtained  for  each  of  4  components  defined  by  the  mode  of 
propagation  of  the  gamma  rays  through  the  barriers: 


Exterior  Barrier 
Propagation 


Interior  Partition 
Propagation 


Designation 


Direct  Beam 
Scattered 
Direct  Beam 
Scattered 


Direct  Beam 
Direct  Beam 
Scattered 
Scattered 


EDB-IDB 

EBS-IDB 

EDB-IBS 

EBS-IBS 


The  components  were  summed  to  give  the  total. 


The  barrier-scattered  components  (i.e.,  EBS-IDB,  EDB-IBS,  and  EBS-IBS) 
were  calculated  by  Monte  Carlo.  The  EDB-IDB  component  was  calculated  by 
analytic  procedures. 

Using  the  Engineering  Method,  protection  factors  were  calculated  for 
the  same  configurations  described  above.  The  Engineering  Method  results  were 
calculated  in  terms  of  a  directr-beam  and  barrier-scattered  exposure  rate 
inside  the  barriers.  This  allowed,  comparison  with  Monte  Carlo  results  on  a 
component  basis. 

Inspection  of  the  Monte  Carlo  results  indicated  that,  as  the  locations 
of  the  interior  partitions  were  varied,  the  Monte  Carlo  protection  factors 
showed  no  consistent  dependence  upon  the  location  of  the  partition.  Overall, 
the  protection  factor  increased  an  average  of  8.0  percent  as  the  radius  in¬ 
creased  from  2.5  to  10.0  ft.  The  location  becomes  more  important  as  the 
barrier  or  partition  thickness  is  increased.  Over  the  range  of  radii  and 
barrier  thicknesses  used  here,  it  was  concluded  that  the  location  of  the 
partition  is  unimportant.  The  neglect  of  the  location  of  the  partition  in 
the  Engineering  Method  is  consistent  with  this  conclusion. 

The  40-psf  partition  increased  the  protection  factors  of  the  various 
barrier  configurations  by  factors  ranging  from  1.97  to  2.47.  The  increases 
were  all  the  same  within  26  percent.  This  indicates  that  some  average 
barrier  reduction  factor  which  is  a  function  only  of  the  partition  thickness 
could  be  used  in  calculating  the  effect  of  the  partition.  This  is  consis¬ 
tent  with  the  Engineering  Method  in  that  the  partition  barrier  reduction 
factor  is  assumed  to  be  a  function  only  of  the  partition  thickness. 

The  Monte  Carlo  and  Engineering  Method  protection  factors  were  com¬ 
pared.  The  average  difference  was  7.7  percent,  with  the  largest  differ¬ 
ence  approximately  22  percent.  The  Engineering  Method  results  were  generally 
higher.  This  same  trend  was  also  noted  in  previous  experimental  studies. 

Average  barrier  reduction  factors  for  the  interior  partitions  were 
calculated  from  the  Monte  Carlo  data  for  each  partition  thickness.  The 
results  were  found  to  be  larger  than  similar  Engineering  Method  data  by 
0,  9,  and  25  percent  for  0-,  20-,  and  40-psf  partitions. 


The  new  barrier  reduction  factors  were  used  in  additional  Engineering 
Method  calculations.  The  results  were  compared  with  Monte  Carlo  protection 
factors  and  found  to  agree  much  better  than  did  the  original  Engineering 
Method  data.  Also,  the  new  results  were  on  the  conservative  side  of  the 
Monte  Carlo  protection  factors. 

It  was  concluded  that  the  assumptions  used  in  the  Engineering  Method 
for  predicting  the  effects  of  interior  partitions  in  fallout  shielding  are 
reasonable.  However,  new  barrier  reduction  factors  should  be  calculated 
for  the  partitions.  At  present,  the  same  barrier  reduction  factor  curves 
are  used  for  both  the  exterior  barriers  and  the  partitions. 

Before  incorporating  new  partition  barrier  reduction  factors,  it  would 
be  necessary  to  extend  the  range  of  the  barrier  and  partition  thicknesses 
and  locations  considered  in  this  study  to  determine  the  exact  range  for 
which  the  assumptions  used  in  the  Engineering  Method  are  valid.  It  would 
also  be  necessary  to  determine  the  exact  dependence  of  on  the  geo¬ 

metric  shape  of  the  partition.  It  should  be  noted  that  these  conclusions 
are  based  on  a  cylindrical  geometry  and  may  not  apply  to  other  geometries. 
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ABSTRACT 


\ 

Protection  factors  were  calculated  by  Monte  Carlo  for  cylindrical 
barriers  with  a  concentric  interior  partition.  The  barrier  thicknesses 
were  20,  AO,  and  80  psf  and  the  partition  thicknesses  were  and  AO  psf. 
The  radius  of  the  barrier  was  10.0  ft  and  the  interior  partition  had  radii 
of  2.5,  5.0, and  10.0  ft.  The  cylindrical  concrete  barriers  were  infinite 
in  height.  The  source  was  an  infinite  plane  60Co  source  which  was  assumed 
to  represent  fallout.  The  Engineering  Method  was  also  used  to  calculate 
protection  factors  for  the  cylindrical  barriers,  and  the  results  were  com¬ 
pared  with  the  Monte  Carlo  data.  ( ' 

i 

The  Monte  Carlo  protection  factors  showed  no  consistent  dependence 
upon  the  radius  of  the  partition.  The  AO-psf  partition  was  found  to  in¬ 
crease  the  protection  factors  of  the  barriers  by  as  much  as  a  factor  of 
2.A7.  The  increases  in  the  protection  factors  due  to  the  AO-psf  partition 
were  all  the  same  within  26  percent  for  the  different  partition  radii  and 
barrier  thicknesses.  The  Monte  Carlo  and  Engineering  Method  protection 
factors  had  an  average  difference  of  7.7  percent  with  the  largest  differ¬ 
ence  being  22  percent.  The  Engineering  Method  results  were  consistently 
higher.  Interior  partition  reduction  factors  were  calculated  from  the 
Monte  Carlo  data.  These  factors  were  found  to  be  as  much  as  25  percent 
higher  than  similar  Engineering  Method  reduction  factors.  Use  of  these 
reduction  factors  in  the  Engineering  Method  gave  protection  factors  which 
agreed  much  better  with  the  Monte  Carlo  protection  factors. 
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SUMMARY 


Protection  factors  were  calculated  for  cylindrical  concrete  barriers 
with  an  interior  partition  using  Monte  Carlo  procedures-  The  results  were 
compared  with  protection  factors  calculated  using  the  Engineering  Method, 
in  order  to  evaluate  the  ability  of  the  Engineering  Method  to  predict  pro¬ 
tection  factors  for  cylindrical  barriers  with  interior  partitions.  The 
Monte  Carlo  results  were  also  analyzed  to  give  a  better  understanding  of 
the  effects  of  the  interior  partition  and  to  aid  in  further  development  of 
the  Engineering  Method. 


The  cylindrical  barriers  were  of  concrete  and  were  infinite  in  height. 
The  radius  of  the  external  barrier  was  10.0  ft,  and  the  interior  partition 
radii  were  2.5,  5.0,  and  10.0  ft.  The  barrier  thicknesses  were  20,  40,  and 
80  psf  (lbs/ft2),  and  the  interior  partition  thicknesses  were  20  and  40  psf. 
The  source  was  an  infinite  plane  60Co  source,  assumed  to  represent  fallout. 


The  results  from  previous  Monte  Carlo  protection  factor  calculations 
for  20-,  40-,  and  80-psf  cylindrical  barriers  without  interior  partitions 
ware  used  in  the  comparisons  to  present  barriers  with  interior  partitions 
of  zero  thickness.  The  barrier  configuration  and  source  were  the  same  as 
those  used  in  the  protection  factor  calculations  with  interior  partitions. 
The  previous  results  were  also  used  to  represent  barriers  with  a  pseudo 
partition  at  a  radius  of  10.0  ft. 


The  flux  energy  and  exposure  rate  angle  distributions  were  also  calcu¬ 
lated  at  a  height  of  3  ft  inside  the  barriers  with  the  interior  partitions. 
These  data  were  obtained  for  each  of  4  components  defined  by  the  mode  of 
propagation  of  the  gamma  rays  through  the  barriers: 


Exterior  Barrier 
Propagation 

Direct  Beam 
Scattered 
Direct  Beam 
Scattered 

The  components  were  s 


Interior  Partition 
Propagation 

Dirp-'t  Beam 
Direct  Beam 
Scattered 
Scattered 

to  give  the  total. 


Designation 

EDB-IDB 

EBS-IDB 

EDB-IBS 

EBS-IBS 
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The  barrier-scattered  components  (i.e.,  EBS-IDB,  EDB-IBS,  and  EBS-IBS) 
were  calculated  by  Monte  Carlo.  The  EDB-IDB  component  was  calculated  by 
analytic  procedures. 

Using  the  Engineering  Method,  protection  factors  were  calculated  for 
the  same  configurations  described  above.  The  Engineering  Method  results  were 
calculated  in  terms  of  a  direct-beam  and  barrier-scattered  exposure  rate 
inside  the  barriers.  This  allowed  comparison  with  Monte  Carlo  results  on  a 
component  basis. 

Inspection  of  the  Monte  Carlo  results  indicated  that,  as  the  locations 
of  the  interior  partitions  were  varied,  the  Monte  Carlo  protection  factors 
showed  no  consistent  dependence  upon  the  location  of  the  partition.  Gverall, 
the  protection  factor  increased  an  average  of  8.0  percent  as  the  radius  in¬ 
creased  from  2.5  to  10.0  ft.  The  location  becomes  more  important  as  the 
barrier  or  partition  thickness  is  increased.  Over  the  range  of  radii  and 
barrier  thicknesses  used  here,  it  was  concluded  that  the  location  of  the 
partition  is  unimportant.  The  neglect  of  the  location  of  the  partition  in 
the  Engineering  Method  is  consistent  with  this  conclusion. 

The  40-psf  partition  increased  the  protection  factors  of  the  various 
barrier  configurations  by  factors  ranging  from  1.97  to  2.47.  The  increases 
were  all  the  same  within  26  percent.  This  indicates  that  some  average 
barrie  •  reduction  factor  which  is  a  function  only  of  the  partition  thickness 
could  be  used  in  calculating  the  effect  of  the  partition.  This  is  consis¬ 
tent  with  the  Engineering  Method  in  that  the  partition  barrier  reduction 
factor  is  assumed  to  be  a  function  only  of  the  partition  thickness. 

The  Monte  Carlo  and  Engineering  Method  protection  factors  were  com¬ 
pared.  The  average  difference  was  7.7  percent,  with  the  largest  differ¬ 
ence  approximately  22  percent.  The  Engineering  Method  results  were  generally 
higher.  This  same  trend  was  also  noted  in  previous  experimental  studies. 

Average  barrier  reduction  factors  for  the  interior  partitions  were 
calculated  from  the  Monte  Carlo  data  for  each  partition  thickness.  The 
results  were  found  to  be  larger  than  similar  Engineering  Method  data  by 
0,  9,  and  25  percent  for  0-,  20-,  and  40-psf  partitions. 
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Tlie  new  barrier  reduction  factors  were  used  in  additional  Engineering 
Method  calculations.  The  results  were  compared  with  Monte  Carlo  protection 
faccors  and  found  to  agree  much  better  than  did  the  original  Engineering 
Method  data.  Also,  the  new  results  were  on  the  conservative  side  of  the 
Monte  Carlo  protection  factors* 

It  was  concluded  that  the  assumptions  used  in  the  Engineering  Method 
for  predicting  the  effects  of  interior  partitions  in  fallout  shielding  are 
reasonable.  However,  new  barrier  reduction  factors  should  be  calculated 
for  the  partitions.  At  present,  the  same  barrier  reduction  factor  curves 
are  used  for  both  the  exterior  barriers  and  the  partitions. 

Before  incorporating  new  partition  barrier  reduction  factors,  it  would 
be  necessary  to  extend  the  range  of  the  barrier  an*  “srtition  thicknesses 
and  locations  considered  in  this  study  to  determine  the  exact  range  for 
which  the  assumptions  used  in  the  Engineering  Method  are  valid.  It  would 
also  be  necessary  to  determine  the  exact  dependence  of  ' (X^)  on  the  geo¬ 
metric  shape  of  the  partition.  It  should  be  noted  that  these  conclusions 
are  based  on  a  cylindrical  geometry  and  may  not  apply  to  other  geometries. 
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I.  INTRODUCTION 


Over  the  past  several  years,  a  small  number  of  experiments  have  been 
performed  to  determine  the  protection  afforded  by  various  structures  with 
interior  partitions  when  exposed  to  a  simulated  fallout  field.1*2*3  The 
results,  expressed  as  protection  factors,  were  compared  with  similar  results 
calculated  with  the  well  known  "Engineering  Method". **  The  Engineering 
Method,  which  was  developed  by  N.  FitzSimons  and  C.  Eisenhauer  from  the 
moments-method  data  reported  in  NBS  Monograph  42, 5  was  found  to  be  in  reason¬ 
able  agreement  with  the  experiments,  although  the  Engineering  Method  pro¬ 
tection  factors  were  generally  high.  However,  it  was  recently  suggested 
that  the  Engineering  Method  for  treating  interior  partitions  is  in  error,® 
and  that  its  success  is  due  to  off-setting  errors  in  the  method. 

The  objective  of  the  present  study  is  to  use  Monte  Carlo  methods  to 
evaluate  and  possibly  develop  further  the  Engineering  Method  treatment  of 
the  interior  partition  problem.  In  particular,  Monte  Carlo  calculations 
should  yield  sufficiently  detailed  information  to  determine  whether  the 
assumptions  and/or  da-ta  used  in  the  Engineering  Method  produce  off-setting 
errors. 

The  geometry  used  for  the  study  consists  of  2  concentric  right  circular 
cylinders  located  normal  to  the  ground  as  shown  in  Figure  1.  The  barriers 
are  assumed  to  be  infinite  in  height;  therefore,  no  roof  contributions  are 
considered.  The  cylindrical  barrier  geometry  was  used  to  be  consistent  with 
that  used  in  previous  Monte  Carlo  studies.7*8  Many  of  the  parameters  used 
in  this  study,  such  as  barrier  thickness,  barrier  radius,  the  energy  and 
angle  groups  for  both  the  source  and  the  results,  and  the  source  term,  were 
the  same  as  those  used  previously. 

This  allowed  the  results  from  the  previous  studies  to  be  used  herein 
for  barriers  with  no  partitions.  The  results  were  also  used  for  barriers 
with  pseudo  partitions.  For  example,  the  results  for  the  80-psf  barrier 
were  used  to  represent  a  40-psf  barrier  with  a  40-psf  partition  whose  outer 
radius  was  equal  to  the  inner  radius  of  the  exterior  barrier. 
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The  barriers  and  partitions  are  composed  of  ordinary  concrete.  The 
source  is  assumed  to  be  an  infinite  plane  of  ^Co(~1.25  MeV)  located  on  the 
ground  surface.  The  detector  was  located  3  ft  above  the  ground  on  the  verti¬ 
cal  axis  of  the  barriers. 

The  parameters  that  were  varied  were  the  exterior  barrier  thickness, 
the  interior  partition  thickness,  and  the  radius  of  the  partition.  The 
values  calculated  for  each  configuration  were  the  flux  energy  distribution 
and  the  exposure  rate  polar  angle  distribution.  The  total  flux  and  exposure 
rate  were  also  determined.  These  data  were  calculated  for  each  of  4  compo¬ 
nents  representing  different  modes  of  propagation  of  the  gamma  rays  through 
the  barrier  and  partition. 

From  these  Monte  Carlo  data,  protection  factors  were  calculated  for 
each  barrier  configuration.  These  factors  were  compared  with  measured  pro¬ 
tection  factors  and  with  protection  factors  calculated  using  the  Engineering 
Method. 

The  Monte  Carlo  calculations  are  described  in  Section  II.  A  descrip¬ 
tion  of  the  source  term  and  its  application  to  the  problem  are  included. 

The  Engineering  Method  for  interior  partition  calculations  is  described  in 
Section  III.  Included  are  the  numerical  values  of  the  various  parameters 
used  in  the  calculations. 

Section  IV  presents  the  Monte  Carlo  results  in  integral  form,  gives 
a  detailed  analysis  of  the  Monte  Carlo  results,  and  compares  them  with  experi¬ 
mental  and  Engineering  Method  results.  The  detailed  energy  and  angle  dis¬ 
tributions  are  given  in  Appendix  A,  Appendix  B  describes  modifications 
made  to  the  GRASS  Monte  Carlo  procedure9  to  make  it  applicable  to  the  inte¬ 
rior  partition  problem  and  to  improve  its  efficiency. 
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II.  MONTE  CARLO  CALCULATIONS 

To  determine  the  protection  factor  of  a  barrier,  it  is  necessary  to 
calculate  1)  the  free-field  exposure  rate  and  2)  the  exposure  rate  inside 
the  barrier.  The  free-field  flux  and  exposure  rate  3  ft  above  an  infinite 
smooth  plane  of  60Co  were  calculated  in  a  previous  study.7  The  results  of 
those  free-field  calculations  were  used  in  this  study. 

The  free-field  flux  was  used  to  develop  a  pseudo  source.  This  source 
was  used  in  the  barrier  penetration  calculations  made  to  determine  the  ex¬ 
posure  rate  inside  the  barriers.  The  following  sub-section  gives  a  brief 
description  of  the  free-field  calculations  described  in  detail  in  Reference 
7.  The  pseudo  source  is  described  in  Section  2.2.  Section  2.3  describes 
the  barrier  penetration  calculations  and  the  protection  factor  calculations. 

2.1  Free-Field  Radiation  Calculations 

The  free-field  calculations  consisted  of  computing  the  flux  energy 
and  angle  distributions  of  both  the  uncollided  and  the  air-  and  ground- 
scattered  gamma  rays  3  ft  above  an  infinite  plane  60CoCVL.25  MeV)  source. 

The  COHORT  Monte  Carlo  procedure10  was  used  for  the  scattered  gamma-ray 
calculations.  The  angle  distribution  of  the  uncollided  gamma  rays  was  com¬ 
puted  analytically  using  the  RRA-65A  procedure.8 

The  polar  angle  distribution  of  the  two  components  (i.e.,  unco'llided 
and  scattered  flux)  were  computed  in  terms  of  the  absolute  number  flux  in 
each  of  18  ten-degree  intervals  (i.e.,  0-10,  10-20,  ...170-180  degrees) 
on  polar  angle.  The  polar  angle  is  measured  with  respect  to  the  normal  to 
the  source  plane.  The  results  were  integrated  over  azimuthal  angle.  The 
scattered  flux  was  sorted  into  8  energy  groups  defined  by  the  limits:  0.04- 
0.06,  0.06-0.10,  0.10-0.18,  0.18-0.30,  0.30-0.50,  0.50-0.75,  0.75-1.00,  and 
1.00-1.25  MeV.  The  free-field  flux  was  converted  to  exposure  rate  using 
the  conversion  factors  from  Reference  11.  A  source  strength  of  1  photon/ 
cm2-sec  was  assumed  for  the  calculations. 
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The  results  of  these  calculations  are  given  in  Table  I  and  consist 
of:  1)  the  flux  energy  and  angle  distribution*  2)  the  total  flux  and  expo¬ 
sure  rate  angle  distribution,  3)  the  total  flux  energy  distribution,  and 
4)  the  total  free-field  flux  and  exposure  rate. 

2.2  Pseudo  Source  Description 

For  the  barrier  penetration  calculations,  it  was  assumed  that  the  energy 
and  angle  distribution  of  the  free-field  flux  is  not  disturbed  by  the  barrier. 
It  has  been  shown  that  the  exposure  rate  inside  a  barrier  is  increased  very 
little  by  gamma  rays  incident  on  the  barrier  more  than  a  few  feet  above  or 
below  the  detector  height.7  Tnus,  it  was  assumed  that  the  flux  incident  at 
any  point  on  the  barrier  could  be  described  by  the  distribution  of  the  free- 
field  flux  at  the  3-ft  detector  height. 

The  energy  and  angle  distribution  of  the  free-field  flux  was  used  to 
generate  a  pseudo  source.  This  source  was  located  on  the  outside  surface  of 
the  exterior  barrier  and  extended  to  a  height  of  23  ft.  The  height  of  the 
barrier  was  46  ft.  The  pseudo  source  is  represented  in  the  barrier  penetra¬ 
tion  calculations  in  the  form  of  probability  distribution  tables  which  are 
defined  in  terms  of  the  energy  and  angle  intervals  given  above  for  the  free- 
field  results.  The  source  is  described  by  a  polar  angle  distribution  for 
each  energy  interval  and  an  energy  probability  distribution.  Due  to  the 
azimuthal  symmetry  of  the  infinite  plane  source,  the  pseudo  source  has  a 
cosine  distribution  on  azimuthal  angle. 

These  probability  distributions  were  used  in  the  Monte  Carlo  calcula¬ 
tions  to  select  the  initial  energy  and  direction  of  each  gamma  ray  started 
from  the  pseudo  source.  The  energy  and  direction  of  the  gamma  rays  were 
assumed  to  be  uniformly  distributed  within  each  angle  and  energy  interval, 
except  the  1.0-1.25  MeV  energy  group.  Inspection  of  Table  I  indicates  that 
approximately  97%  of  the  flux  in  this  group  is  due  to  uncollided  gamma  rays 
that  have  an  energy  of  1.25  MeV.  For  this  reason,  approximately  97%  of  the 
gamma  rays  within  this  energy  group  were  given  an  initial  energy  of  1.25 
MeV.  The  remainder  were  evenly  distributed  over  the  energy  interval. 
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2.3  Barrier  Penetration  Calculations 


The  quantities  calculated  in  the  barrier  penetration  calculations 
include:  1)  the  total  flux,  2)  the  total  exposure  rate,  3)  the  flux  energy 
distribution,  and  4)  the  exposure  rate  polar  angle  distribution.  The  f.lux 
energy  distribution  was  calculated  for  eight  energy  groups  and  the  exposure 
rate  angle  distributions  for  18  ten-degree  polar  angle  groups.  The  energy 
and  angle  groups  are  the  same  as  those  used  for  the  free-field  calculations.. 
Each  of  these  quantities  was  calculated  for  a  point  detector  located  at  a 
height  of  3  ft  on  the  centerline  of  the  concentric  cylindrical  concrete  bar¬ 
riers.  The  results  were  integrated  over  azimuthal  angle. 

The  calculations  were  made  for  the  12  possible  combinations  of  the  fol¬ 
lowing  parameters: 

1)  Partition  thicknesses  of  20  and  40  psf 

2)  Exterior  barrier  ithicknesses  of  20,  40,  and  80  psf 

3)  Partition  radii  of  2.5  and  5.0  ft. 

The  exterior  barrier  had  an  outer  radius  of  10  ft. 


Each  of  the  quantities  described'  above  was  calculated  separately  for 
each  of  4  components  defined  by  the  mode  of  gamma-ray  propagation  through 
the  barriers: 


Exterior  Barrier 
Propagation 


Partition 

Propagation 


Designation 


Direct  Beam 
Scattered 
Direct  Beam 
Scattered 


Direct  Beam 
Direct  Beam 
Scattered 
Scattered 


EDB-IDB 

EBS-IDB 

EDB-IBS 

EBS-IBS 


The  four  components  were  summed  to  give  the  total  of  each  quantity  calculated. 
Knowing  the  total  exposure  rate  inside  each  barrier  and  the  free-field  expo¬ 
sure  rate,  the  protection  factor  for  each  of  the  barriers  may  be  determined. 
The  protection  factor  is  defined  as  the  free-field  exposure  rate  divided  by 
the  exposure  rate  inside  the  barrier. 
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The  barrier  penetration  calculations  were  performed  by  separate  methods 
for  the  barrier-scattered  components  (i.e.,  EBS-IDB,  EDB-IBS,  and  EBS-IBS) 
and  the  direct-beam  component  (EDB-IDB) .  The  calculations  are  described 
separately  below.  The  results  of  the  penetration  calculations  and  the  pro¬ 
tection  factors  are  given  in  Section  IV. 

2.3.1  Barrier-Scattered  Components 

The  barrier-scattered  components  were  computed  using  a  modified  version 
of  the  GRASS  Monte  Carlo  procedure.9  The  modifications  were  made  to  include: 
1)  a  generalized  geometry  capable  of  treating  basements  as  well  as  interior 
partitions,  2)  the  separation  of  the  exposure  rate  inside  the  barrier  into 
components  that  are  defined  by  the  mode  of  propagation  of  the  gamma  rays 
through  the  barrier,  and  3)  biasing  techniques  to  improve  the  statistical 
accuracy  of  the  EDB-IBS  and  EBS-IBS  components. 

The  modified  version  of  GRASS  is  designated  GRAS1.  A  description  of 
the  modifications  is  given  in  Appendix  B  and  a  complete  description  of  GRASS 
is  given  in  Reference  9.  GRAS1  will  be  briefly  described  here  as  it  pertains 
to  the  interior  partition  problem. 

GRAS1  calculates  the  barrier-scattered  flux  energy  and  exposure  rate 
polar  angle  distribution  at  a  point  inside  concentric  cylindrical  barriers. 
The  gamma-ray  source  is  a  pseudo  source  located  on  the  outside  surface  of 
the  exterior  barrier.  This  pseudo  source  is  defined  by  probability  distri¬ 
butions  obtained  from  the  energy  and  angle  distributions  of  the  incident 
flux,  and  by  the  total  number  of  photons  incident  upon  the  barrier.  (See 
Section  2.2.) 

Once  the  initial  gamma-ray  energy,  direction,  and  spatial  position 
are  determined  from  the  proper  probability  distributions,  Monte  Carlo  tech¬ 
niques  are  used  to  trace  the  path  of  the  gamma  ray  as  it  scatters  within 
the  defined  geometry.  Each  gamma  ray  is  traced  until  it  reaches  a  specified 
minimum  energy  or  weight,  has  a  specified  number  of  collisions,  or  escapes 
from  the  defined  geometry.  A  statistical  estimate  of  the  scattered  flux 
reaching  a  point  detector  is  made  from  each  collision  point.  The  sum  of 


these  estimates  divided  by  the  number  of  gamma  rays  traced  gives  an  estimate 
of  the  scattered  flux  per  photon  incident  upon  the  barrier.  The  total  scat¬ 
tered  flux  is  obtained  by  multiplying  this  quantity  by  the  number  of  incident 
photons. 

The  flux  is  recorded  for  each  of  the  three  barrier-scattered  components 
in  terms  of  eight  energy  groups  and  the  exposure  rate  polar  angle  distribu¬ 
tion  in  18  ten-degree  intervals.  All  results  are  integrated  with  respect 
to  azimuthal  angle.  Also  calculated  are  the  total  flux  and  exposure  rate. 

To  determine  if  GRAS1  was  working  properly*  GRAS1  calculations  were 
made  of  the  energy  and  angle  distributions  inside  a  40-psf  exterior  cylindri¬ 
cal  barrier  with  an  adjacent  40-psf  partition.  The  results  were  compared 
with  the  results  from  the  previous  Monte  Carlo  calculations  for  an  80-psf 
barrier.  The  total  exposure  rate  was  found  to  agree  within  2.5%  while  the 
total  flux  agreed  within  1.5%.  A  comparison  of  the  exposure  rate  polar  angle 
distribution  is  shown  in  Figure  2.  The  results  are  seen  to  agree  quite  well, 
but  the  GRAS1  results  (based  on  only  1000  histories)  are  not  as  smooth  as 
the  GRASS  Monte  Carlo  data,  which  are  for  a  much  larger  number  of  histories. 
Similar  agreement  was  found  for  the  flux  energy  distribution.  Comparisons 
made  for  barriers  having  a  total  thickness  of  20  and  40  psf  gave  similar 
results.  It  was  concluded  that  GRAS1  was  working  properly. 

From  these  comparisons,  it  was  also  concluded  that  the  results  of  the 
previous  study  could  be  used  as  a  basis  of  comparison  in  the  present  study. 
Any  differences  found  between  the  results  for  the  barrier  with  and  without 
partitions  are  caused  by  the  partition  and  not  by  differences  in  the  methods 
of  calculation. 

GRAS1  calculations  were  made  in  steps  of  500  histories  per  problem  to 
give  some  indication  of  the  statistical  accuracy  of  the  calculations.  To 
get  from  16,000  to  18,000  collisions  for  each  configuration,  it  was  neces¬ 
sary  to  run  two,  three,  and  four  problems  for  the  80-,  40-,  and  20-psf 
barrier  configurations,  respectively.  The  number  of  collisions  per  history 
decreases  as  the  barrier  thickness  decreases.  The  results  from  like  pro¬ 
blems  were  averaged  to  give  the  final  results.  The  standard  deviation  of 


EXPOSURE 


1 


the  total  barrier-scattered  exposure  rate  for  the  20-psf  barriers  was  found 
to  range  from  9  to  20  percent  with  the  average  standard  deviation  being 
approximately  14  percent. 

The  accuracy  of  the  results  for  the  40t  and  80-psf  barriers  should  be 
about  the  same,  because  approximately  the  same  number  of  collisions  was 
used  for  each  configuration. 

All  cross  sections  used  in  the  calculations  were  taken  from  Reference 
12.  The  exposure  rate  conversion  factors  were  the  same  as  those  used  in  the 
free-field  radiation  calculations. 

2.3.2  Direct-Beam  Component 

The  direct-beam  component  (EDB-IDB)  was  calculated  using  existing 
analytical  procedures  RRA-65A®  and  RRA-67.7  The  calculations  were  done  in 
2  parts  which  are  defined  by  the  mode  of  propagation  of  the  gamma  rays  before 
reaching  the  barrier:  1)  gamma  rays  incident  on  the  barrier  after  having 
at  least  one  collision  in  the  air  and/or  ground,  and  2)  gamma  rays  incident 
on  the  barrier  before  undergoing  a  collision . 

The  RRA-67  procedure  was  used  for  the  first  part.  This  procedure  com¬ 
putes  attenuation  factors  for  the  barrier  and  partition  for  each  energy  and 
angle  group  of  the  air-  and  ground-scattered  free-field  flux  given  in  Table 
I.  It  multiplies  these  factors  by  the  corresponding  flux  values  given  in 
the  table.  The  procedure  then  performs  a  #iux-to-exposure  rate  conversion 
and  sums  over  energy  and  angle  to  obtain  the  total  exposure  rate. 

The  second  part  of  the  EDB-IDB  component  is  one  of  the  most  important 
components  of  the  total  exposure  rate.  It  was  calculated  by  the  RRA-65A 
procedure  which  integrates  over  the  entire  infinite  plane  source  and  applies 
exponential  attenuation  to  the  flux  from  each  differential  source  area. 

The  results  were  then  converted  to  exposure  rate  and  added  to  the  results 
of  part  1  of  the  calculation  to  obtain  the  total  direct-beam  component. 
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III.  ENGINEERING  METHOD  CALCULATIONS 

The  Engineering  Method  for  predicting  the  protection  factor  of  a 
structure  is  described  in  detail  in  Reference  4  (OCD-PM- 100-1) .  The  brief 
description  given  in  this  section  describes  the  method  for  handling  problems 
involving  interior  partitions.  No  roof  contributions  are  considered.  The 
curves  and  charts  given  in  OCD-PM-100-1  are  based  on  moments-method  data  for 
1.12-hr  fallout  given  in  NBS-42.5  Corresponding  moments-method  data  for  60Co 
are  given  in  Appendix  B  of  NBS-42.  The  necessary  curves  for  60Co  were  gen¬ 
erated  from  these  data  for  use  in  the  Engineering  Method  protection  factor 
calculations.  The  geometry  factors  used  for  60Co  were  those  calculated  by 
LeDoux  from  the  NBS-42  moments-method  data  and  reported  by  Kaplan  in  Ref¬ 
erence  2. 


Protection  factors  were  calculated  using  the  Engineering  Method  for 
each  of  the  barrier  configurations  described  in  Section  II.  The  equations 
used  for  the  calculations  are  described  below.  The  sources  of  data  and  the 
numerical  values  of  the  many  parameters  required  in  the  calculations  are 
also  given. 


3.1  Description  of  Method 

The  Engineering  Method  for  predicting  the  protection  factor  of  a  cylin¬ 
drical  barrier  of  infinite  height  with  an  interior  partition  consists  of 
1)  determining  a  reduction  factor  assuming  no  partition,  and  2)  multiplying 
this  reduction  factor  by  an  interior  partition  barrier  reduction  factor. 

The  protection  factor  is  the  reciprocal  of  this  product. 

The  Engineering  Method  assumes  that: 

1.  The  protection  afforded  by  the  partition  is  independent  of  the 
solid  angle  subtended  by  the  partition. 

2.  The  partition  barrier  reduction  factor  is  independent  of  the 
shape  of  the  partition. 

3.  The  same  barrier  reduction  factor  curves  may  be  used  for  both 
the  exterior  barrier  and  the  partition. 

The  partition  barrier  reduction  factor  is  a  function  only  of  the  partition 
thickness. 
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The  Engineering  Method  equation  for  the  reduction  factor  for  a  cylin¬ 
drical  barrier  with  no  partition  is 


*  Be(Xe*H){fGd(V  +Gu(V][1-Sw(Xe)]} 


+  Be(Xe,H){tGs(Ujl)  +  V«u)HSw.(Xe)][E]} 


where  Be(Xe>H)  is  the  barrier  reduction  factor  at  a  detector  height  H  for  a 

vertical  slab  of  thickness  X  , 

e’ 

Gd is  the  directional  response  function  for  direct  radiation, 

W  is  the  directional  response  function  for  air-scattered  radiation, 

Gs^J^  ls  the  directional  response  function  for  barrier-scattered 
radiation  from  below  the  detector  plane, 

is  ihe  directional  response  function  for  barrier-scattered 
radiation  from  above  the  detector  plane, 

Sw^Xe^  ls  a  wei8hting  factor  that  approximates  the  fraction  of  the  radia¬ 
tion  reaching  the  detector  which  is  scattered  in  the  barrier, 

and  E  is  the  shape  factor  for  barrier-scattered  radiation  (E  =  tt/2  for  a 
cylinder) . * 3 


The  directional  response  functions  are  a  function  of  the  solid  angle 
fraction,  w,  which  is  defined  as  the  solid  angle  subtended  by  the  exterior 
barrier  at  the  detector,  divided  by  2tt.  The  subscripts  u  and  £  are  used  to 
denote  the  solid  angle  fractions  above  and  below  the  detector,  respectively. 
For  a  cylinder,  o>  =  1.0-cosineG  where  6  is  the  angle  formed  between  the 
centerline  and  a  line  connecting  the  top  or  bottom  of  the  barrier  to  the 
detector. 

The  protection  factor  for  a  barrier  with  a  partition  is  given  by  the 
equation 


P  -  1.0 

f  [Rf][Bi(xi)] 

where  B^)  is  the  barrier  reduction  factor  for  a  partition  of  thickness  X  . 
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Since  the  weighting  factor,  S  (X  •) ,  approximates  the  fraction  of  radia¬ 
tion  reaching  the  detector  which  has  been  scattered  in  the  barrier,13  the 
first  term  on  the  right  hand  side  of  the  reduction  factor  equation  may  be 
interpreted  as  the  direct-beam  component  and  the  second,  term  as  the  barrier- 
scattered  component.  To  allow  direct  comparisons  with  the  Monte  Carlo 
results,  the  Engineering  Method  results  are  presented  in  terms  of  the  direct- 
beam  and  barrier-scattered  components  in  addition  to  the  total  protection 
factor. 

3.2  Parameter.  Data  and  Results 

The  figures  and  charts  given  in  OCD-PM-lOO-1  are  based  on  the  moments- 
method  data  given  in  Reference  4  (NBS-42)  for  1.12-hr  fallout.  Corresponding 
moments-method  data  for  60Co  are  given  in  Appendix  B  of  NBS-42.  The  necessary 
curves  for  60Co  were  generated  from  these  data.  The  moments-method  data  used 
to  evaluate  the  parameters  in  the  Engineering  Method  calculations  are: 

Parameter  Figure  Reference 

Be,  B±  Fig.  B25  Ref.  5 

G,,  G  ,  G  Fig.  2  Ref.  2 

Sw  Fig.  B23  Ref.  5 

It  is  often  very  difficult  to  duplicate  protection  factor  calculations 
for  even  a  simple  structure  because  of  minor  differences  in  reading  the 
numerical  value  of  the  various  parameters  from  the  curves.  Thus,  to  docu¬ 
ment  fully  the  calculations  reported  here,  the  numerical  values  of  each  of 
the  parameters  used  are  given  in  Table  II. 

The  results  of  the  Engineering  Method  protection  factor  calculations 
are  given  in  Table  III.  The  direct-beam  and  barrier-scattered  components 
are  expressed  in  units  of  R/hr  per  source  photon/cm2-sec.  The  results  are 
based  on  a  free-field  exposure  rate  of  6.055  x  10"6  R/hr  obtained  from  the 
moments-method  data  in  NBS  Monograph  42. 
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Table  II.  Numerical  Value  of  Parameters  for 
Engineering  Method  Calculations 


Parameter 

X(psf) 

Parameter 

to 

u 

20 

AO 

80 

0.026 

0.713 

B 

0.625 

0.380 

0.1A3 

G 

0.496 

0.278 

e 

s 

B. 

0.625 

0.380 

G 

0.088 

1 

u 

S 

0.327 

0.505 

0.699 

G  . 

0.642 

w 

d 

Table  III.  Engineering  Method  Data  for  Concentric  Cylindrical 
Barriers  Exposed  to  Infinite  Plane  60Co  Source 


(R/hr  per  source  photon/cm2-sec) 


Radiation 

Barrier  Thickness 

Component 

20 

40 

.  80 

P-psf 

Interior  Partition 

Direct  Beam 

1.859-6* 

8.324-7 

1.901-7 

Scattered 

1.504-6 

1.413-6 

7.363-7 

Total 

3.363-6 

2.245-6 

9.264-7 

Protection** 

Factor 

1.800 

2.696 

6.535 

20-psf 

Interior  Partition 

Direct  Beam 

1.162-6 

5.201-7 

1.188-7 

Scattered 

9.405-7 

8.833-7 

4.605-7 

Total 

2.103-6 

1.404-6 

5.793-7 

Protection 

Factor 

2.879 

4.314 

10.450 

40-psf 

Interior  Partition 

Direct  Beam 

7.070-7 

3.162-7 

7.228-8 

Scattered 

5.719-7 

5.369-7 

2.800-7 

Total 

1.279-6 

8.531-7 

3.523-7 

Protection 

Factor 

4.734 

7.097 

17.190 

*  Read  1.859-6  as  1.859  x  l(f6  . 

**  Dq  =  6.055  x  10~6  (R/hr  per  source  photon/cm2~sec) 
based  on  moments  method  data  from  Reference  5. 
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IV.  RESULTS  AND  ANALYSIS 

The  results  of  the  Monte  Carlo  barrier  penetration  calculations  are 
summarized  in  Table  IV.  The  exposure  rate  is  given  for  each  of  the  four 
components  described  in  Section  2.3.  The  detector  is  at  a  height  of  3  ft 
on  the  centerline  of  the  concentric  cylindrical  barriers.  The  barriers 
have  exterior  thicknesses  of  20,  40,  and  80  psf  and  partition  thicknesses 
of  20  and  40  psf.  The  results  are  given  for  partition  radii  of  both  2.5 
and  5.0  ft.  the  outer  radius  of  the  exterior  barrier  is  10.0  ft,  and  the 
infinite  plane  60Co  source  on  the  ground  surface  has  a  strength  of  1  photon/ 
cm2-sec. 

Also  given  in  Table  IV  is  the  protection  factor  of  each  barrier.  The 
protection  factor  is  defined  as  the  free-field  exposure  rate  divided  by  the 
total  exposure  rate  inside  the  barrier.  The  total  exposure  rate  was  obtained 
by  summing  the  results  of  the  four  components. 

Given  in  Appendix  A  are  the  flux  energy  distributions  and  the  exposure 
rate  angle  distributions  in  differential  form  for  each  component.  The  re¬ 
sults  are  given  for  each  of  8  energy  groups  and  18  ten-degree  polar  angle 
groups.  The  results  are  integrated  over  azimuthal  angle. 

In  a  previous  study,7  protection  factors  were  calculated  for  20-,  40-, 
and  80-psf  cylindrical  barriers  using  the  COHORT  Monte  Carlo  procedure.10 
The  barriers  had  no  partitions,  but  the  barrier  configuration,  source,  and 
detector  were  the  same  as  used  in  this  study.  These  results  are  used  as  a 
basis  for  comparison  in  determining  the  effects  of  the  partitions  and  also 
as  illustrations  of  limiting  conditions  for  barriers  with  partitions,  as 
described  below.  These  results  are  given  in  Table  V  for  20-,  40-,  and  80-psf 
barriers  and  as  interpolated  for  a  60-psf  barrier.  The  exposure  rates  are 
given  for  the  direct-beam  component  and  the  barrier-scattered  component. 

The  scattered  component  corresponds  to  the  sum  of  the  scattered  components 
in  Table  IV.  The  protection  factors  are  given  for  each  barrier. 


Table  IV,  Penetration  Data  for  Concentric  Cylindrical  Barriers 
Exposed  to  Infinite  Plane  60Co  Source 
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Read  1.294-6  as  1.294  x  10-6  . 

D  =  5.989-6  (R/hr  per  source  photon/cm2-sec)  . 
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TABLE  V.  Penetration  Data  for  Simple  Cylindrical  Barriers  Exposed 
to  Infinite  Plane  60 Co  Source 


(E/hr  per  source  photon/cm2-sec) 


Radiation 

Component 

- 

Barrier  Thickness  (psf) 

20 

40 

60* 

~  80 

Direct  Beam 

2.332-6** 

1.294-6 

7.135-7 

4.053-7 

Scattered 

7.343-7 

8.214-7 

7.370-7 

4.848-7 

Total 

3.066-6 

2.116-6 

1.451-6 

8.901-7 

Protection 

Factor 

1.954 

2.830 

4,125 

6.734 

*  The  scattered  component  for  60  psf  was  obtained  by  interpolation  of  the 
20-,  40-,  and  80-psf  data. 

**  Read  2.332-6  as  2.332  x  10"6 . 


The  40-,  60- ,  and  80-psf  data  are  used  in  the  analysis  as  barriers 
with  pseudo  partitions  having  a  radius  of  approximately  10  ft.  This  is 
possible  by  considering  the  80-psf  barrier,  with  the  10-ft  outer  radius,  to 
consist  of  two  adjacent  40-psf,  concentric  barriers.  Therefore,  the  80-psf 
data  is  the  same  as  that  for  a  40-psf  barrier  with  a  40-psf  partition. 
Similarly,  the  60-psf  data  are  used  for  20-  and  40-psf  barriers  with  40- 
and  20-psf  partitions,  respectively. 

4.1  Analysis  of  Components 

Before  considering  the  protection  factors,  it  is  advantageous  to 
examine  first  the  effects  of  the  partitions  on  the  various  exposure  rate 
components  considered  in  the  protection  factor  calculations.  It  is  also 
of  interest  to  determine  the  effects  of  the  partitions  on  .the  energy  and 
angle  distributions  of  the  flux  inside  the  barriers. 

It  is  suggested  in  Reference  6  that  the  Engineering  Method  over¬ 
estimates  the  direct-beam  component  and  underestimates  the  scattered 
components  by  approximately  equal  amounts.  Thus,  it  was  concluded  that, 
although  the  overall  results  were  reasonably  good,  the  Engineering  Method 
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appeared  Co  be  .in  error.  Comparisons  of  these  components  of  the  Monte  Carlo 
and  Engineering  Method  data  presented  in  Tables  IV  and  V  and  Table  III, 
respectively,  indicate  just  the  opposite.  As  an  example,  for  a  40-psf  bar¬ 
rier  with  no  partition,  the  Engineering  Method  overpredicts  the  importance 
of  the  scattered  components  by  as  much  as  62  percent  while  underpredicting 
the  direct  beam  by  as  much  as  40  percent. 

In  the  Engineering  Method,  the  importance  of  the  barrier-scattered  com¬ 
ponent  is  assumed  to  be  independent  of  the  partition  thickness.  However, 
the  Monte  Carlo  data  (Table  IV)  indicate  that  this  component  increases  from 
24,  39,  and  55  percent  of  the  total  exposure  for  a  20-,  40-,  and  80-psf 
barrier,  respectively,  -to  53,  61,  and  70  percent  for  barriers  with  40-psf 
partitions . 

For  the  40-psf  partition,  the  Engineering  Method  was  found  to  under- 
predlct  the  Monte  Carlo  data  by  as  much  as  18  percent  for  the  20-psf  barrier. 
The  largest  overprediction  was  only  21.5  percent.  Despite  the  improvement 
in  the  comparisons  for  the  barriers  with  partitions,  the  large  differences 
noted  for  the  barriers  without  partitions  are  unexplained  and  should  be  in¬ 
vestigated.  From  the  results  of  this  and  previous  studies,7*8  no  conclusive 
explanation  for  these  differences  can  be  made. 

In  the  analysis  of  the  overall  results,  it  is  important  to  have  some 
indication  of  the  importance  of  the  various  barrier-scattered  components. 
First,  consider  the  barriers  with  the  20-psf  partitions.  The  EBS-IDB  compo¬ 
nent  ranges  from  12  to  32  percent  of  the  total  exposure  rate  for  the  20-  to 
80-psf  barriers  with  both  the  2.5-  and  5.0-ft  partition  radius.  For  both 
partition  radii,  the  EDB-IBS  components  range  from  15  to  17,  18  to  21,  and 
10  to  12  percent  for  the  20-,  40-,  and  80-psf  barriers,  respectively.  For 
the  2.5-ft  partition  radius,  the  EBS-IBS  component  accounted  for  13  and  21 
percent  of  the  total  exposure  inside  the  20-  and  80-psf  barrier,  respectively. 
Corresponding  percentages  of  only  8  and  17  percent  were  found  with  a  parti¬ 
tion  radius  of  5.0  ft.  This  component  shows  a  definite  decrease  in  importance 
as  the  partition  radius  is  increased. 
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For  the  40-psf  interior  partition,  the  EDB-IBS  and  EBS-IBS  components 
show  a  definite  increase  in  their  importance  to  the  total  exposure.  This 
is  as  would  be  expected  and,  in  some  cases,  the  increase  is  as  large  as  a 
factor  of  2.  However,  the  overall  differences  among  the  various  components 
remains  about  the  same  as  for  the  20-psf  partition. 

For  both  the  20-  and  40-psf  partition,  no  individual  barrier-scattered 
component  was  consistently  an  extremely  important  component.  The  importance 
of  a  scattered  component  may  vary  greatly  from  configuration  to  configuration. 

Comparisons  were  made  of  the  flux  energy  distributions  and  the  exposure 
rate  polar  angle  distributions  of  the  sum  of  all  radiation  components  inside 
the  cylindrical  barriers  both  with  and  without  partitions.  Figure  3  shows 
the  angle  distribution  inside  a  20-psf  barrier  both  without  a  partition  and 
with  a  40-psf  partition  at  2.5-ft  radius.  The  angle  distribution  is  not  as 
peaked  around  90  degrees  for  the  configuration  with  the  partition  as  it  is 
without  the  partition.  The  difference,  which  is  slight,  was  found  for  both 
20-  and  40-psf  partitions  and  all  barrier  thicknesses.  No  appreciable  dif¬ 
ference  was  found  between  the  angle  distributions  inside  the  barriers  with 
a  5-ft  radius  partition  and  that  inside  corresponding  barriers  without  par¬ 
titions. 

The  energy  spectra  inside  a  40-psf  barrier  both  with  and  without  a 
partition  are  shown  in  Figure  4.  The  40-psf  partition  has  a  radius  of  5.0 
ft.  The  energy  spectrum  is  slightly  harder  for  the  configuration  with  the 
partition.  This  was  also  found  true  for  the  other  configurations  considered 
(see  Appendix  A) . 

4.2  Effect  on  Protection  Factor  of  Partition  Location 

Given  in  Tables  IV  and  V  are  the  protection  factors  for  20-  and  40-psf 
barriers  with  20-  and  40-psf  partitions  having  radii  of  2.5,  5.0,  and  10.0 
ft.  These  protection  factors  are  plotted  versus  partition  radius  in  Figure 
5.  Inspection  of  the  protection  factors  indicates  no  consistent  dependence 
upon  the  partition  radius,  with  the  exception  that  the  protection  factors 
are  always  higher  for  the  10-ft  radius  than  for  the  2.5-ft  radius.  For  the 
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Fig.  3.  Angle  Distribution  of  the  Scattered  Exposure  Rate  Inside  a  20-psf  Barrier  with  and 
without  40-psf  Partition 
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Fig.  4.  Scattered  Flux  Energy  Distribution  Inside  a  40-psf  Barrier  with  and  without  40-psf  Partition 
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80-psf  barrier,  the  protection  factors  were  determined  only  for  the  2.5-  and 
5.0-ft  partition  radius.  These  data  indicate  that  the  protection  factors 
for  the  5.0-ft  radius  are  higher  than  those  for  the  2.5-ft  radius  by  11  and 
14  percent  for  the  20-  and  40rpsf  partitions,  respectively.  The  Engineering 
Method  is  independent  of  the  partition  location. 

Examination  of  the  barrier-scattered  components  indicated  that  only  the 
EBS-IBS  component  shows  any  consistent  dependence  upon  the  partition  location. 
This  component  varies  from  8  to  31  percent  of  the  total  exposure  rate.  With 
one  exception,  the  exposure  rates  due  to  the  EBS-IBS  component  for  the  2.5-ft 
partition  radius  are  30  to  40  percent  higher  than  for  the  5.0-ft  configuration. 
For  the  20-psf  barrier  with  the  40-psf  partition,  the  exposure  rate  is  13 
percent  less  for  the  2.5-ft  data.  The  2.5-ft  data  appears  to  be  somewhat  low. 
The  data  for  the  10.0-ft  configuration  were  calculated  for  the  total  barrier- 
scattered  component  only. 

It  is  concluded  that  for  barrier  thicknesses  and  configurations  where 
the  EBS-IBS  component  Is  the  principle  component,  the  protection  factor  is 
reduced  as  the  partition  is  placed  closer  to  the  detector  position.  However, 
for  barrier  and  partition  thicknesses  and  radii  within  the  range  used  in  this 
study,  the  location  of  the  partition  can  be  neglected.  This  conclusion  is 
in  agreement  with  the  assumptions  used  in  the  Engineering  Method.  Any  error 
introduced  into  the  results  by  neglecting  the  partition  location  is  on  the 
conservative  side. 

4.3  Effect  on  Protection  Factor  of  Barrier  and  Partition  Thickness 

The  Monte  Carlo  protection  factors  are  plotted  versus  exterior  barrier 
thickness  in  Figure  6.  Also  shown  are  experimental  results  which  will  be 
discussed  in  Section  4.4. 

With  no  partition,  the  protection  factor  for  an  80-psf  barrier  is  a 
factor  of  3.-45  larger  than  that  for  a  20-psf  barrier.  With  40-psf  partitions, 
the  protection  factor  for  the  80-psf  barrier  is  as  much  as  a  factor  of  4.15 
greater  than  for  the  20-psf  barrier.  This  increase  of  approximately  20  per¬ 
cent  appears  to  be  due  mainly  to  the  changes  in  the  EDB-IBS  component.  As 
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Fig.  6.  Protection  Foctor  vs.  Exterior  Barrier  Thickness  for  Concentric  Cylindrical  Barriers 


indicated  in  Table  IV,  this  component  increases  approximately  27.5  percent 
more  for  a  20-psf  barrier  than  for  an  80-psf  barrier  when  the  partition 
thickness  is  increased  from  20  to  40  psf .  The  same  effect  was  found  for  par¬ 
tition  radii  of  2.5  and  5.0  ft.. 

With  40-psf  partitions,  the  protection  factors  for  the  various  barrier 
thicknesses  and  partition  radii  ranged  from  a  factor  of  1.97  to  4.27  larger 
than  the  protection  factors  for  corresponding  barriers  with  no  partitions. 
Overall,  these  factors  show  a  slight  increase  for  the  40-psf  partition  thick¬ 
ness  as  the  exterior  barrier  thickness  is  increased.  However,  as  the  barrier 
thickness  was  increased  from  20-  to  80-psf,  there  was  less  than  a  26  percent 
difference  between  the  largest  and  smallest  increase  in  the  protection  factors 
due  to  the  addition  of  the  partition.  This  small  variation  suggests  that 
perhaps  some  average  barrier  reduction  factor  could  be  used  for  the  partitions 
It  would  be  a  function  only  of  the  partition  thickness,  a  fact  which  is  con¬ 
sistent  with  the  neglect  of  the  exterior  barrier  thickness  in  the  Engineering 
Method. 

4.4  Comparisons  of  Monte  Carlo,  Experimental,  and  Engineering 
Method  Protection  Factors 

It  is  assumed  in  the  Engineering  Method  that  the  protection  factor  for 
a  barrier  with  a  partition  is  independent  of  location  of  the  partition.  It 
is  also  assumed  that  the  barrier  reduction  factor  for  the  partition  is  a 
function  only  of  the  partition  thickness.  Using  Monte  Carlo  data  in  the 
preceding  sections,  these  assumptions  were  found  to  be  sound.  However, 
verification  of  the  assumptions  gives  no  indication  of  the  overall  ability 
of  the  Engineering  Method,  to  predict  the  protection  factor  of  a  barrier  with 
a  partition.  Comparisons  between  the  Monte  Carlo  and  Engineering  Method 
protection  factors  will  be  used  to  determine  the  ability  of  the  Engineering 
Method  to  predict  protection  factors.  The  Monte  Carlo  protection  factors 
are  compared  with  similar  measured  data  before  using  them  as  a  basis  of 
comparison. 

Technical  Operations,  Inc.,  conducted  an  experiment  in  which  scale 
models  were  exposed  to  a  60Co  source.2  One  of  the  scale  models  was  a  cylin¬ 
drical  steel  barrier  2  ft  in  height  and  1  ft  in  radius  with'  an  interior 
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partition,  0.5  ft  in  radius.  The  °°Co  source  simulated  a  uniformly  contami¬ 
nated  infinite  plane  of  gamma  radiation  and  was  comparable  to  the  source  term 
assumed  in  the  present  penetration  calculations.  Because  the  barriers  were 
not  of  infinite  height,  the  exposure  rates  inside  the  barriers  contained  a 
skyshine  component  which  is  not  present  in  the  Monte  Carlo  results.  To  allow 
a  direct  comparison,  the  skyshine,  through  the  open  top,  was  removed  by  sub¬ 
tracting  an  Engineering  Method  estimate  of  this  component. 

The  adjusted  measured  results  for  a  detector  height  of  3  in.  were  used 
to  calculate  protection  factors  for  model  barriers  of  various  thicknesses. 

The  3-in.  data  were  used  because  the  solid  angle  fractions  subtended  at  this 
height  were  closest  to  those  for  which  the  calculations  were  made.  The  re¬ 
sulting  protection  factors  are  compared  with  the  Monte  Carlo  protection 
factors  in  Figure  6. 

The  results  agree  within  3.5  percent  for  the  barriers  with  no  interior 
partitions.  For  the  barriers  with  a  20-psf  partition,  the  Monte  Carlo  results 
for  the  5.0-ft-radius  partition  agree  with  the  measured  results  wither.  20 
percent  and  the  Monte  Carlo  results  for  the  2.5-ft-radius  partition  agree 
within  10  percent.  The  5.0-ft-radius  partition  was  expected  to  give  the  best 
comparison  since  the  two  geometries  are  more  comparable.  The  total  attenua¬ 
tion  afforded  by  the  partitions  was  found  to  be  generally  15  to  20  percent 
less  for  the  measured  data  than  for  the  Monte  Carlo  data.  Differences  of 
this  magnitude  could  be  due  to  the  different  barrier  configurations  and 
sources  used. 

Figures  7,  8,  and  9  show  comparisons  of  the  protection  factors  at  a 
height  of  3  ft  inside  the  20-,  40-,  and  80-psf  barriers,  respectively,  with 
0— ,  20-,  and  40-psf  partitions  calculated  by  Monte  Carlo  and  by  the  Engi¬ 
neering  Method.  Included  in  the  comparisons  are  the  results  of  the  Engi¬ 
neering  Method  calculations  using  a  barrier  reduction  factor  based  on  the 
Monte  Carlo  data.  These  special  calculations  will  be  described  below. 

The  agreement  between  the  Monte  Carlo  and  the  regular  Engineering 
Method  calculations  is  reasonably  good  in  both  shape  and  magnitude  with  the 
largest  difference  being  approximately  22  percent  for  the  40-psf  barrier 
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Fig.  7.  Protection  Factor  vs.  Partition  Thickness  for  Concentric  Cylindrical  Barriers  (Partition 
Radius  =2.5  ft ) 
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Fig.  8.  Protection  Foctor  vs.  Partition  Thickness  for  Concentric  Cylindrical  Barriers  (Partition 
Radius  =  5.0  ft ) 


31 


PARTITION  THICKNESS  (psf) 


Fig.  9.  Proteclion  Factor  vs.  Partition  Thickness  for  Concentric  Cylindrical  Barriers  (Partition 
Radius  =  10.0  ft ) 
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with  a  40-psf  partition  at  a  radius  of  5.0  ft.  The  results  for  the  20-psf 
barrier  with  a  20-psf  partition  agree  within  1  percent.  The  average  differ¬ 
ence  for  all  barrier  configurations  is  7.7  percent.  Overall,  the  agreement 
is  generally  better  for  the  barriers  with  the  thinner  interior  partitions. 
However,  as  pointed  out  in  Section  4.1,  the  best  agreement  between  the  .Monte 
Carlo  and  the  Engineering  Method  -results  for  the  barrier-scattered  component 
was  found  for  the  thicker  interior  partitions. 

The  -increased  disagreement  observed  as  the  interior  barrier  becomes 
more  important  indicates:  1)  that  the  barrier  factors  assumed  for  the  interior 
partitions  are  incorrect,  or  2)  that  the  method  to  account  for  the  partition 
is  incorrect.  As  indicated  by  all  previous  experimental  studies  1»2*3  as  well 
as  in  Figures  7,  8,  and  9,  the  protection  factors  predicted  by  the  Engineering 
Method  are  nonconservative.  Thus,  it  seems  desirable  to  modify  the  Engineering 
Method  to  give  improved  agreement  or, at  least,  a  conservative  prediction. 

Since  the  preceding  analysis  suggests  that  the  assumptions  used  in  the 
Engineering  Method  are  not  unreasonable,  it  was  decided  to  try  to  improve 
the  interior  partition  barrier  factors.  The  Monte  Carlo  reduction  factors 
for  a  given  barrier  with  a  partition  was  divided  by  that  for  the  same  barrier 
with  no  partition  to  obtain  Monte  Carlo  barrier  reduction  factors.  As  a  re¬ 
sult,  six  new  barrier  reduction  factors  were  determined  for  each  partition 
thickness  (i.e.,  one  for  each  of  the  three  external  barrier  thicknesses  and 
for  both  the  2.5  and  5.0-ft  radius  partition). 

As  suggested  in  Section  4.3,  an  average  barrier  reduction  factor  was 
determined  for  each  partition  thickness  by  averaging  the  six  new  barrier 
reduction  factors.  The  average  Monte  Carlo  barrier  reduction  factor,  B_. '  (X^) , 
was  found  to  be  larger  than  the  original  factors,  B^(X^),  (see  Section  III) 
by  0,  9,  and  25  percent  for  0-,  20-,  and  40-psf  partitions,  respectively. 

It  should  be  noted,  however,  that  the  shape  of  the  partition  is  accounted 
for  in  the  terms.  The  shape  of  the  partition  is  not  considered  in 

the  Engineering  Method. 

Figures  7,  8,  and  9  include  a  comparison  with  Engineering  Method 
results  obtained  by  replacing  B^(X^)  with  The  agreement  is  somewhat 
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improved  between  the  Monte  Carlo,  and  the  new  Engineering  Method  results,  and, 
in  addition,  the  new  Engineering  Method  results  are  generally  conservative. 

Recall  that  the  Engineering  Method  accounts  for  the  partition  by  adjusting 
the  reduction  factor  for  a  barrier  with  no  partition  by  a  barrier  reduction 
factor  which  is  a  function  only  of  partition  thickness.  The  same  reduction 
factor  curves  are  used  for  both  the  barrier  and  the  partition.  The  results 
and  improvements  observed  above  indicate  that  the  Engineering  Method  may  be 
adequate  if  separate  data  are  used  for  the  interior  partition  reduction  factors. 

Changing  the  Engineering  Method  to  incorporate  separate  barrier  reduction 
factor  data  for  the  partitions  would  involve  only  the  addition  of  B^'(X^)  data 
to  an  existing  barrier  reduction  factor  chart  (i.e..  Chart  1,  Reference  4). 
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V.  CONCLUSIONS 

The  main  objectives  of  this  study  were  to  determine  the  effects  of 
partitions  in  fallout  shielding  and  the  ability  of  the  Engineering  Method 
to  predict  such  effects.  Monte  Carlo  procedures  were  used  to  calculate 
the  exposure  rate  inside  20- ,  40- ,  and  80-psf  cylindrical  barriers  with 
0- ,  20-,  and  40-psf  interior  partitions.  Tne  radius  of  the  cylindrical 
barriers  was  10.0  ft.  The  barriers  were  exposed  to  an  infinite  plane  60Co 
source.  Engineering  Method  calculations  were  also  made  to  determine  the 
effects  of  the  interior  partitions. 

The  final  results  of  the  calculations  were  expressed  in  the  form  of  a 
protection  factor  for  each  barrier.  Ah  analysis  of  the  Monte  Carlo  data 
was  made  to  determine  the  effects  of  the  partition.  The  Monte  Carlo  and 
Engineering  Method  protection  factors  were  compared  to  determine  the  abi¬ 
lity  of  the  Engineering  Method  to  predict  such  effects. 

The  Monte  Carlo  protection  factors  showed  no  consistent  dependence 
upon  the  location  of  the-  partition-.  Overall,  the  protection  factor  in¬ 
creased  an  average  of  8.0  percent  as  the  radius  increased  from  2.5  to  10.0 
ft.  The  location  becomes  more  important  as  the  barrier  or  partition  thick¬ 
ness  is  increased.  Over  the  range  of  radii  and  barrier  thicknesses  used 
here,  it  was  concluded  that  the  location  of  the  partition  is  unimportant. 

The  neglect  of  the  location  of  the  partition  in  the  Engineering  Method  is 
consistent  with  this  conclusion. 

The  40-psf  partition  increased  the  protection  factors  of  the  various 
barrier  configurations  by  factors  ranging  from  1.97  to  2.47.  The  increases 
were  all  the  same  within  26  percent.  This  indicates  that  some  average 
barrier  reduction  factor  which  is  a  function  only  of  the  partition  thick¬ 
ness  could  be  used  in  calculating  the  effect  of  the  partition.  This  is 
consistent  with  che  Engineering  Method  in  that  the  partition  barrier  re¬ 
duction  factor  is  assumed  to  be  a  function  only  of  the  partition  thickness. 

The  Monte  Carlo  and  Engineering  Method  protection  factors  were  compared. 
The  average  difference  was  1.1  percent,  with  the  largest  difference  approxi¬ 
mately  22  percent.  The  Engineering  Method  results  were  generally  higher. 

This  same  trend  was  also  noted  in  previous  experimental  studies. 


Average  barrier  reduction  factors  for  the  interior  partitions  were 
calculated  from  the  Monte  Carlo  data  for  each  partition  thickness.  The 
results  were  found  to  be  larger  than  similar  Engineering  Method  data  by 
0,  9,  and  25  percent  for  0-,  20-,  and  40-psf  partitions. 

The  new  barrier  reduction  factors  were  used  in  additional  Engineering 
Method  calculations.  The  results  were  compared  with  Monte  Carlo  protection 
factors  and  found  to  agree  much  better  than  did  the  original  Engineerng 
Method  data.  Also,  the  new  results  were  oh  the  conservative  side  of  the 
Monte  Carlo  protection  factors. 

It  was  concluded  that  the  assumptions  used  in  the  Engineering  Method 
for  predicting  the  effects  of  interior  partitions  in  fallout  shielding  are 
reasonable.  However,  new  barrier  reduction  factors  should  be  calculated 
for  the  partitions.  At  present,  the  same  barrier  reduction  factor  curves 
are  used  for  both  the  exterior  barriers  and  the  partitions. 

Before  incorporating  new  partition  barrier  reduction  factors,  it  would 
be  necessary  to  extend  the  range  of  the  barrier  and  partition  thicknesses 
and  locations  considered  in  this  study  to  determine  the  exact  range  for 
which  the  assumptions  used  in  the  Engineering  Method  are  valid.  It  would 
also  be  necessary  to  determine  the  exact  dependence  of  ’  (X_^)  on  the  geo¬ 
metric  shape  of  the  partition.  It  should  be  noted  that  these  conclusions 
are  based  only  on  cylindrical  geometries  and  may  not  apply  to  other  geometries. 
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APPENDIX  A 


'Gamnia-Ray  Energy  and  Angle  Distributions 
Inside  Concentric  Cylindrical  Barriers  Exposed  to 
Infinite  Plane  ^°Co  Source 

Table  Page 

A1  Flux  Energy  Distribution  Inside  20- ,  40- ,  38 

and'  80-psf,  Barriers  with  a  20-psf  Partition 
at  a  2.5-ft  Radius 

A2.  Flux  Energy  Distribution  Inside  20-,  40-,  39 

and  80-psf  Barriers  with  a  20-psf  Partition 
at  a  5.0-ft  Radius 

A3  Flux  Energy  Distribution  Inside  20- ,  40- ,  40 

and  80-psf  Barriers  with  a  40-psf  Partition 
at  a  2.5-ft  Radius 

A4  Flux  Energy  Distribution  Inside  20- ,  40- ,  41 

and  80-psf  Barriers  with  a  40-psf  Partition 
at  a  5.0-ft  Radius 

A5  Exposure  Rate  Angle  Distribution  Inside  a  42 

20-psf  Barrier  with  a  20-psf  Partition 

A6  Exposure  Rate  Angle  Distribution  Inside  a  43 

40-psf  Barrier  with  a  20-psf  Partition 

A7  Exposure  Rate  Angle  Distribution  Inside  a  44 

80-psf  Barrier  with  a  20-psf  Partition 

A8  Exposure  Rate  Angle  Distribution  Inside  a  45 

20-psf  Barrier  with  a  40-psf  Partition 

A9  Exposure  Rate  Angle  Distribution  Inside  a  46 

40-psf  Barrier  with  a  40-psf  Partition 

A10  Exposure  Rate  Angle  Distribution  Inside  a  47 

80-psf  Barrier  with  a  40-psf  Partition 


The  distributions  are  listed  in  terms  of  the  upper 
bounds  of  the  energy  and  angle  groups;  the  lower 
bound  is  0.04  MeV  for  the  first  energy  group  and  0° 
for  the  first  angle  group. 


NOTE: 
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TABLE  Al.  Flux  Energy  Distribution  Inside  20-,  40-,  and  80-psf 
Barriers  with  a  20-psf  Partition  at  a  2. 5-ft  Radius 


’  2  2 
.  '(photons/cm  -sec  per, source  photon/cm  -sec) 


Energy 

(Mev) 

Radiation  Component 

Total 

EDB-IBS 

EBS-IDB 

EBS-IBS 

EDB-IDB 

20-psf  Barrier  Thickness 

0.06 

4.857-03* 

'2.346-04 

2.946-02 

0.000-00 

,  3.455-02 

0.10 

2.589-02 

6.888-03 

1.429-01 

0.000-00 

1.757-01 

0;18 

9.412-02 

2.657-02 

2.089-01 

0.000-00 

3.296-01 

0.30 

9.505-02 

3.228-02 

1.034-01 

0.000-00 

2.307-01 

0.50 

3.893-02 

3.092-02 

1.070-01 

0.000-00 

1.769-01 

0.75 

6.724-02 

4.110-02 

4.947-02 

0.000-00 

1.578-01 

1.00 

1.936-02 

4.651-02 

1.099-02 

0.000-00 

7.686-02 

1.25 

5.144-02 

4.035-02 

0.000-00 

5.953-01 

6.871-01 

Total 

3.969-01 

2.248-01 

6.525-01 

5.953-01 

1.869-00 

40-psf  Barrier  Thickness 

0.06 

4.296-03 

3.232-04 

1.306-02 

■0.000-00 

1.768-02 

0.10 

5.636-02 

7.191-03 

1.065-01 

0.000-00 

1.701-01 

0.18 

9.762-02 

2.590-02 

1.487-01 

0.000-00 

2.722-01 

0.30 

1.018-01 

2.957-02 

9.015-02 

0.000-00 

2.215-01 

0.50 

4.513-02 

3.235-02 

6.569-02 

0.000-00 

1.432-01 

0.75 

2.9  92-02 

4.795-02 

2.460-02 

0.000-00 

1.025-01 

1.00 

3.131-02 

4.225-02 

3.960-02 

0.000-00 

1.132-01 

1.25 

5.494-02 

4.027-02 

0.000-00 

3.233-01 

4.185-01 

Total 

4.215-01 

2.258-01 

4.885-01 

3.233-01 

1.459-00 

80-psf  Barrier  Thickness 

0.06 

7.089-04 

1.624-04 

6.121-03 

0.000-00 

6.992-03 

0.10 

8.620-03 

3.370-03 

8.517-02 

0.000-00 

9.716-02 

0.18 

1.142-02 

1.558-02 

1.098-01 

0.000-00 

1.368-01 

0.30 

1.629-02 

1.334-02 

5.027-02 

0.000-00 

7.990-02 

0.50 

9.387-03 

1.515-02 

2.369-02 

0.000-00 

4.823-02 

0.75 

1.343-02 

3.000-02 

3.551-02 

0.000-00 

7.894-02 

1.00 

1.575-02 

2.681-02 

2.219-03 

0.000-00 

4.478-02 

1.25 

0.000-00 

4.264-02 

1.957-04 

9.911-02 

1.419-01 

Total 

7.563-02 

1.470-01 

3.131-01 

9.911-02 

6.347-01 

*  Read  4.857-03  as  4.857  x  10"3. 
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TABLE  A2.  Flux  Energy  Distribution  Inside  20-,  4Q-,  and  80-psf 
Barriers  «Lth  a  20-psf  Partition  at  a  5.0-  ft  Radius 

2  2 
(photons/cm  -sec  per  source  photon/cm  -sec) 


Energy 


Radiation  Component  ; 


(Mev) 


EDB-IBS  EBS-IDB  EBS-IBS  EDB-IDB  Total 


20-psf  Barrier  Thickness 


0.06 

1.792-02 

0.10 

9.246-02 

0.18 

1.273-01 

0.30 

8.717-02 

0.50 

6.682-02 

0.75 

5.617-02 

1.00 

4.436-02 

1.25 

3.407-02 

Total 

5.263-01 

1.918-04  2.469-02 
6.781-03  1.044-01 
2.784-02  1.466-01 
2.831-02  6.976-02 
3.136-02  6.047-02 
4.419-02  2.760-02 
3.875u02  3.592-03 
4.414-02  0. 000-00 

2.215-01  4.373-01 


0.000-00  4.282-02 
0.000-00  2.036-01 
0.000-00  3.019-01 
0.000-00  1.852-01 
0.000-00  1.586-01 
0.000-00  1.279-01 
0.000-00  8.671-02 
5.953-01  6.735-01 

5.953-01  1.780-00 


40-nsf  Barrier  Thickness 


0.06 

3.919-03 

2.408-04 

0.10 

3.804-02 

8.313-03 

0.18 

4.817-02 

2.577-02 

0.30 

4.081-02 

3.012-02 

0.50 

2.827-02 

3.619-02 

0.75 

4.026-02 

5.088-02 

1.00 

2.676-02 

4.121-02 

1.25 

4.994-02 

4.697-02 

Total 

2.762-01 

2.396-01 

2.005-02 

0.000-00 

2.421-02 

1.145-01 

0.000-00 

1.609-01 

1.659-01 

0.000-00 

2.398-01 

7.164-02 

0.000-00 

1.426-01 

5.551-02 

0.000-00 

1.200-01 

1.965-02 

0.000-00 

1.108-01 

6.503-03 

0.000-00 

7.447-02 

0.000-00 

3.233-01 

4.202-01 

4.539-01 

3.233-01 

1.293-00 

80-psf  Barrier  Thickness 


0.06 

9.319-04 

1.472-04 

0.10 

6.712-03 

3.950-03 

0.18 

1.355-02 

1.361-02 

0.30 

1.521-02 

1.205-02 

0.50 

6.828-03 

1.628-02 

0.75 

1.043-02 

3.376-02 

1.00 

4.507-03 

2.918-02 

1.25 

1.249-02 

2.606-02 

Total 

7.068-02 

1.350-01 

3.746-03 

0.000-00 

’4.825-03 

3.613-02 

0.000-00 

4.679-02 

5.74  7-02 

0.000-00 

8.463-02 

4.600-02 

0.000-00 

7.326-02 

2.087-02 

0.000-00 

4.398-02 

2.796-02 

0.000-00 

7.215-02 

9.910-07 

0.000-00 

3.369-02 

0.000-00 

9.911-02 

1.376-01 

1.921-01 

9.911-02 

4.969-01 
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TABLE  A3.  Flux  Energy  Distribution  Inside  20-,  40-,  and  '80-psf 
Barriers  with  a  40-psf  Partition  at  a  2.5-  ft  Radius 


2  2 
(photons/cm  -sec  per  source  photon/cm  -sec) 


Energy 

(Mev) 

Radiation  Component' 

Total 

EDB-IBS 

EBS-IDB 

EBS-IBS 

EDB-IDB 

20-psf  Barrier  Thickness 

0.06 

5.801-03 

7.660-06 

.8.566-03 

0.000-00 

1.437-02 

0.10 

8.806t02 

8.600-04 

6.256-02 

0.000-00 

1.515-01 

0.18 

1.489-01 

4.948-03 

9.367-02 

0.000-00 

2.475-01 

0.30 

8,434-02 

6.939-03. 

5.666-02 

.0.000-00 

1.479-01 

0.50 

7.389-02 

9.415-03 

4.009-02 

0.000-00 

1.234-01 

0.75 

1.305-01 

1.728-02 

2.742-02 

0.000-00 

1.752-01 

1.00 

7.656-02 

2.164-02 

4.032-04 

0.000-00 

9.860-02 

1.25 

5.077-02 

2.852-02 

0.000-00 

3.233-01 

4.026-01 

Total 

6.589-01 

8.961-02 

2.893-01 

3.233-01 

1.361-00 

40-psf  Barrier  Thickness 

0.06 

3.546-03 

9.400-06 

7.144-03 

0.000-00 

1.070-02 

0.10 

3. 576-02 

7.923-04 

5.096-02 

0.000-00 

8.751-02 

0.18 

4.794-02 

4.711-03 

1.332-01 

0.000-00 

1.859-01 

0.30 

4.-383-02 

8.992-03 

8.066-02 

0.000-00 

1.335-01 

0.50 

3.587-02 

1.: 169-02 

5.981-02 

0.000-00 

1.074-01 

0.75 

4,179-02 

1.990-02 

5.159-02 

0.000-00 

1.133-01 

1.00 

4.202-02 

2.524-02 

2.968-02 

0.000-00 

9.694-02 

1.25 

3.198-02 

2.459-02 

2.690-04 

1.783-01 

2.351-01 

Total 

2.827-01 

9.594-02 

4.134-01 

1.783-01 

9.703-01 

80-psf  Barrier  Thickness 

■ 

0.06 

1.399-03 

4.485-06 

2.308-03 

0.000-00 

3.711-03 

0.10 

9.232-03 

5.465-04 

4.464-02 

0.000-00 

5.442-02 

0.18 

1.322-02 

3.316-03 

7.315-02 

0.000-00 

8.969-02 

0.30 

2.116-02 

4.379-03 

5.392-02 

0.000-00 

7.946-02 

0.50 

6.579-03 

5.693-03 

1.573-02 

0.000-00 

2.800-02 

0.75 

1.460-02 

8.892-03 

5.481-03 

0.000-00 

2.897-02 

1.00 

4.783-03 

1.399-02 

3.609-02 

0,000-00 

5.486-02 

1.25 

1.779-02 

1.488-02 

0.000-00 

5.530-02 

8.797-02 

Total 

8.877-02 

5.171-02 

2.312-01 

5.530-02 

4.270-01 
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TABLE  A4.  Flux  Energy  Distribution  Inside  20-,  40-,  and  80-psf 
Barriers  with  a  40-psf  Partition  at  a  5.0r  ft  Radius 


2  2 
(photons/cm. -sec  per  source  photon/cm  -sec) 


Energy 

(Mev), 

Radiation  Component- 

Total 

EDB-IBS 

EBS-IDB 

EBS-IBS- 

EDB-IDB 

20-psf  Barrier  Thickness 

0.06 

7.758-03 

6.521-06 

7.360-03 

0.000-00 

1.512-02 

0.10 

6.046-02 

7.236-04 

6.672-02 

0.000-00 

1.2.79-01 

0.18 

1.273-01 

4.402-03 

9.939-02 

0.000-00 

2.311-01 

0.30 

1.083-01 

7.044-03 

4.885-02 

0.000-00 

1.642-01 

0.50 

6.842-02 

1.053-02 

3.377-02 

0.000-00 

1.127-01 

0.75 

8.142-02 

1.907-02 

2.250-02 

0.000-00 

1.230-01 

1.00 

6.932-02 

2.091-02 

8.606-03 

0.000-00 

9.884-02 

1.25 

6.803-02 

2.617-02 

8.642-03 

3.233-01 

4.261-01 

Total 

5.912-01 

8.889-02 

2.958-01 

3.233-01 

1.299-00 

40-psf  Barrier  Thickness 

0.06 

6.824-03 

7.637-06 

8.875-03 

0.000-00 

1.571-02 

0.10 

3.801-02 

8.768-04 

7.677-02 

0.000-00 

1.157-01 

0.18 

5.116-02 

5.120-03 

1.058-01 

0.000-00 

1.621-01 

0.30 

4.283-02 

7.437-03 

6.366-02 

0.000-00 

1.139-01 

0.50 

4.167-02 

1.091-02 

4.810-02 

0.000-00 

1.007-01 

0.75 

3.858-02 

1.957-02 

4.044-02 

0.000-00 

9.859-02 

1.00 

3.131-02 

2.024-02 

2.441-02 

0.000-00 

9.596-02 

1.25 

4.732-02 

3.896-02 

5.530-03 

1.783-01 

2.701-01 

Total 

3.176-01 

1.031-01 

3.736-01 

1.783-01 

9.727-01 

80-psf  Barrier  Thickness 

0.06 

1.800-03 

3.327-06 

6.197-03 

0.000-00 

.8.000-03 

0.10 

9.793-03 

4.620-04 

3.181-02 

0.000-00 

4.207-02 

0.18 

1.616-02 

2.560-03 

4.394-02 

0.000-00 

6.266-02 

0.30 

1.600-02 

3.921-03 

2.458-02 

0.000-00 

4.450-02 

0.50 

6.833-03 

7.405-03 

1.626-02 

0.000-00 

3.050-02 

0.75 

1.382-02 

1.214-02 

8.697-03 

0.000-00 

3.466-02 

1.00 

1.382-02 

1.162-02 

2.221-02 

0.000-00 

4.765-02 

1.25 

5.465-03 

1.500-02 

0.000-00 

5.530-02 

7.577-02 

Total 

8.372-02 

5.313-02 

1.536-01 

5.530-02 

3.458-01 

TABLE  A5.  Exposure  Rate  Angle  Distribution  Inside 
20-psf  Barrier  With  a  20-psf  Partition 


i 
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TABLE  A6.  Exposure  Rate  Angle  Distribution  Inside 
40-psf  Barrier  with  a  20-psf  Partition 
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TABLE  A7.  Exposure  Rate  Angle  Distribution  Inside 
80-psf  Barrier  with  a  20-psf  Partition 
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APPENDIX  B 

GRAS1  Monte  Carlo  Procedure 


The  GRAS1  Monte  Carlo  Procedure  is  a  modified  version  of  the  GRASS 
Monte  Carlo  Procedure. 9  GRASS  was  written  in  FORTRAN  for  an  IBM  1130 
computer  to  calculate  the  barrier-scattered  flux  energy  and  exposure  angle 
distribution  at  a  point  inside  a  cylindrical  barrier  exposed  to  gamma  rays 
from  a  plane  source.  GRASS  was  written  specifically  for  fallout  shielding 
calculations.  A  complete  description  of  the  procedure  along  with  utiliza¬ 
tion  instructions  is  given  in  Reference  9.,  The  GRAS1  procedure  is  basically 
the  same  as  GRASS;  therefore,  only  the  changes  will  be  described  here. 

The  modifications  incorporated  in  GRAS1  include  1)  a  generalized 
cylindrical  geometry  capability,  2)  the  separation  of  the  results  into 
components,  and  3)  biasing  techniques  to  improve  the  statistical  accuracy 
of  the  components.  These  changes  were  made  to  make  the  GRASS  procedure 
applicable  to  fallout  shielding  calculations  where  interior  partitions  and/ 
or  basements  are  considered. 

GRASS,  has  a  fixed  geometry  consisting  of  a  single  cylindrical  barrier 
located  upright  above  a  j)lane  gamma  source.  The  GRAS1  procedure  will  accept 
any  geometry  made  up  of  material  regions  which  are  defined  by  two  types  of 
boundaries:  1)  concentric  cylinders  and  2)  planes  perpendicular  to  the 
centerline  of  the  cylindrical  boundaries .  The  centerline  is  perpendicular 
to  the  source  plane.  GRAS1  is  limited  to  eight  material  regions  and  eight 
boundaries.  The  materials  are  defined  by  two  sets  of  cross  section  data 
and  up  to  five  material  densities.  The  mathematical  relationships  describing 
the  boundaries  and  regions  were  taken  from  the  COHORT  procedure.10 

In  fallout  shielding  calculations  involving  interior  partitions,  it 
is  instructive  to  classify  the  radiation  arriving  at  the  detector  according 
to  the  regions  in  which  scattering  occurred.  In  GRAS1,  the  results  are 
divided  into  3  components  defined  by  the  mode  of  propagation  of  the  gamma 
rays  through  the  barrier  and  partition: 
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1.  Gamma  rays  having  the  initial  collision  in  a  partition 
(EDB-IBS) . 

2.  Gamma  rays  having  collisions  only,  iii  the  exterior  barrier 
(EBS-IDB) . 

3.  Gamma  rays  having  collisions  in  both  the  exterior  barrier 
and  the  partition  (EBS-IBS). 


These  results  are  summed  to  give  the  total  barrier-scattered  flux 
and  exposure  rate. 

The  number  of  gamma-ray  collisions  occurring  in  a  partition  that  is 
located  within  a  somewhat  thicker  barrier  may  be  very  small  compared  to  the 
number  of  collisions  occurring  in  the  barrier.  If  similar  statistical 
accuracy  is  to  be  obtained  for  each  component,  approximately  the  same  number 
of  collisions  should  occur  in  both  the  barrier  and  the  partition.  Therefore, 
special  biasing  was  incorporated  to  improve  the  statistical  accuracy  of  the 
EDB-IBS  and  EBS-IBS  components.  A  different  type  of  biasing  was  used  for 
each  component. 

To  improve  the  statistical  accuracy  of  the  EDB-IBS  component,  it  is 
necessary  to  increase  the  number  of  gamma  rays  passing  directly  through 
the  exterior  barrier.  This  is  done  by  using  the  exponential  transformation 
to  increase  the  mean  free  path  of  the  gamma  rays  before  the  first  collision. 
The  exponential  transformation  is  well  known  and  has  frequently  been  used 
in  shielding  calculations  involving  shields  several  mean-free  paths  in 
thickness. 

-In  brief,  a  pseudo  total  macroscopic  cross  section  is  computed  using 
the  equation 


EX(E)  =  ET(E)[1-ET] 


where  E^,(E)  is  the  total  macroscopic  cross  section  for  the  exterior  barrier 
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at  energy  E,  and  ET  is  an  input  parameter  <1.0.  This  pseudo  cross,  section 
is  then  used  instead  of  Z^,(E)  to  compute  the  random  path  length  to  the  first 
collision. 

The  biasing  introduced  in  the  gamma-ray  weight  is  removed  by  adjusting 
the  weight  of  the  gamma  ray.  If  the  gamma  ray  has  its  initial  collision 
in  the  exterior  barrier,  the  new  gamma-ray  weight  is  given  by  the  equation 

v  (vS  _ 

w.  =  W-I _  e-D[ET(E)-ET(E)3 

ZT(  E) 


where  W  is  the  original  weight  and  D  is  the  distance  to  the  collision  point. 

If  the  gamma  ray  passes  through  the  exterior  barrier  before  having  a  collision, 
the  weight  is  given  by  the  equation 


.  ,  ^e-St£T(E)-ET(E)] 


where  S  is  the  path  length  of  the.  gamma  ray  through  the  exterior  barrier. 

The  exponential  transformation  is  applied  to  a  fraction  of  the  source 
gamma  rays  which  are  selected  at  random.  This  fraction  is  an  input  para¬ 
meter.  Biasing  is  applied  only  to  gamma  rays  that  would  intersect  an 
inside  partition  if  the  random  path  length  selected  were  large  enough. 

To  improve  the  statistical  accuracy  of  the  EBS-IBS  component,  the 
number  of  gamma  rays  that  have  a  collision  in  the  barrier  and  go  to  the 
partitiqn  before  undergoing  the  next  collision  must  be  increased.  This 
is  done  by  splitting  a  fraction  of  the  gamma  rays  which  have  collisions 
in  the  barrier  into  two  gamma  rays.  This  fraction  is  an  input  parameter 
in  the  program,  and  the  gamma  rays  to  be  split  are  selected  at  random. 

Trial  scattering  angles  for  one  of  the  two  gamma  rays  are  selected  from 
the  Klein-Nishina  formula  until  an  angle  is  selected  such  that  the  gamma 
ray  projected'  through  this  angle  will  intersect  the  partition.  Knowing 
the  gamma  ray  will  hit  the  partition,  it  is  forced  to  scatter  at  the  selected 
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angle,  directly  to  the  partition.  At  the  position  on  the  partition  surface 
where  the  gamma  enters,  a  sub-history  is  formed  and  the  gamma  ray  is  traced 
until  it  is  killed.  Sub-histories  are  killed  by  the  same  procedures  as.  the 
original  histories. 


The  weight  of  the  sub-history  is  calculated  from  the  equation 


W* 


where  W  is  the  weight  of  the  original  particle, 

N  is  the  number  of  scattering  angles  selected  before  a  successful 
angle  was  calculated,  and 

p  is  the  number  of  mean  free  paths  between  the  collision  point  and 
the  inside  partition. 


The  original  gamma- ray  history  continues  at  the  end  of  the  sub-r 
history.  Its  new  weight  is 
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